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We generate entangled states of an ensemble of 5 x 10* *^Rb atoms by optical quantum nonde- 
molition measurement. The resonator-enhanced measurement leaves the atomic ensemble, prepared 
in a superposition of hyperfine clock levels, in a squeezed spin state. By comparing the resulting re- 
duction of quantum projection noise (up to 8.8(8) dB) with the concomitant reduction of coherence, 
we demonstrate a clock input state with spectroscopic sensitivity 3.0(8) dB beyond the standard 
quantum limit. 
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Atomic clocks [H-ll] and atom interferometers 4] are 
reaching the standard quantum limit (SQL) of precision 
0, H, [J, set by the quantum projection noise inherent 
in measurements on a collection of uncorrelated parti- 
cles. In the canonical Ramsey interferometer with A^o 
particles, a quantum mechanical phase is converted into 
occupation probabilities for two states and read out as 
a population difference A^ between them. Entanglement 
can reduce the projection noise AN by redistributing it 
to another variable that does not directly affect the ex- 
periment precision. The resulting "squeezed spin state" 
|7h16| can be used as an input state to an interferometer 
to overcome the SQL 0, i, fl Q . 

Formally, the system can be described by an ensemble 
spin vector S = ^ that is the sum over the (pseudo-) 
spins Si of the individual (spin-1/2) particles [5|-l7|. The 
ensemble spin S with (S"^) = S{S + 1) can take on val- 
ues in the range < 5 < S'o, where 5o = No/2. For a 
given 5*, the minimum variance AS'^ of Sz = N/2 for an 
unentangled state is realized by the coherent spin state 
(CSS), and is given by A^^^gg = 5/2 = |(S)|/2, where 
it is assumed that the mean ensemble spin vector (S) lies 
in the xy-plane. A spin state can be defined as squeezed 
if it satisfies Cc = SAS'^/KS)! < 1 (entanglement cri- 
terion 0,[ill), or Cm = 2AS^SiJ\{S)\^ < 1 (criterion 
for metrological gain Q , where Sin is the initial spin 
of the uncorrelated ensemble before the squeezing). 
represents the increase in the squared signal-to-noise ra- 
tio |(S)|^/A5^ over the value 2Sin for the initial uncor- 
related state. Since |(S)| < Sin, we have Cc < Cm, i-e. 
metrological gain guarantees entanglement. 

The process utilized for spin squeezing can reduce |(S)| 
below the initial spin Sin before the squeezing, thereby re- 
ducing the minimum variance AS*^ that is consistent with 



an unentangled state Therefore, measurements of 



both spin noise ASz and average spin length after squeez- 
ing I (S) I are necessary to verify spin squeezing or quantify 
metrological gain. While reduction of spin noise alone has 
sometimes been referred to as "spin squeezing" [l^, El] 
or "number squeezing" H, 2^, we take spin squeezing 
to require at least demonstrated entanglement, Co < 1, 
although we are primarily interested in metrological gain. 

Cm < 1. 



(a) 

Probe 
Trap 



Cavity 



Photod elector 



Filter 



-500 lis 




(b) s^RbDJine 
52n 



|F=2,m=0> 



5^5, 



3.57 GHz 
']6.8GHz 

i |F=1,m=0) 



FIG. 1: (a) Experimental setup, (b) Atomic level structure, 
(c) Experimental sequence. Timing of probe pulses (solid 
line) and microwave pulses (dashed line) in preparation and 
readout of a squeezed state, tt designates a composite n pulse 
[2H . Various procedures are inserted at X, as described in 
the text, to measure the CSS variance, measure the noise of 
a spin component other than Sz, or operate a clock. A-C 
illustrate semiclassical probability distribution functions for 
the Gaussian states discussed in the text. 



Spin noise has been modified by atomic collisions 
|l9, 20, 22] and by absorption of squeezed light In 
dilute atomic systems, quantum nondemolition (QND) 



measurements with 




have re- 



duced the projection noise of rotating |17| and station- 
ary [1^ spins. Spin squeezing has been achieved with 
two ions (8|, and spectroscopic sensitivity further im- 
proved with a maximally entangled state of three ions 
,9]. Recently, spin squeezing with a Bose-Einstein con- 
densate (BEC) in a multiple-well potential has been re- 
ported 24]. Demonstrated metrological gains over the 
SQL include Cm^ — 3.2(1) dB in the three-ion system [§]; 
Cm^ ^4 dB by light- induced squeezing within individual 
atoms of larg e spin s — 3 [lH; and Cm^ = 3.8(4) dB for 
the BEC [23. 

In this Letter, we demonstrate the generation of 
squeezed spin states of 5 x 10'* trapped ^''Rb atoms on 
an atomic-clock transition by resonator-aided QND mea- 
surement with a far-detuned light field, as proposed by 
Kuzmich, Bigelow, and Alandel [13]. We verify the en- 
tanglement by comparing the observed reduction in pro- 
jection noise below that of a coherent spin state (up to 



2 



dB) with the accompanying reduction in clock sig- 
nal, and achieve a 3.0(8) dB improvement in precision 
over the SQL. 

The light-induced spin squeezing presented here re- 
quires strong ensemble-light coupling 

[13, EMI (large 

collective cooperativity [2l|). This is achieved by means 
of a near-confocal optical resonator with, at the 27: /k = 
780 nm wavelength of the probe light, a finesse J- — 
5.6(2) X fO^, a linewidth k ^ 2tt x f.01(3)MHz, and 
a mode waist w = 56.9(4) ym at the atoms' position, 
corresponding to a maximal sing le-atom cooperativity 
770 = 24J'/(7rfc2iy2) ^ 0.203(7) [IJ. Our experiments are 
performed on an ensemble containing up to N^, = 5 x 10"* 
laser-cooled ^''Rb atoms optically trapped inside the res- 
onator in a standing wave of 851-nm light (Fig. [IJ. 

One resonator mode is tuned 3.57(1) GHz to the blue 
of the |52S'i/2,-F = 2) ^ |52P3/2,F' = 3) transition in 
^''Rb, such that the atomic index of refraction results 
in a mode frequency shift lo that is proportional to the 
population difference N = N2 — Ni between the hy- 
perfine clock states |1) = |5^S'i/2, = 1, mj? = O) and 
|2) = \5'^Si/2,F — 2,mF — O). The transmission of a 
probe laser tuned to the slope of this mode thus directly 
measures Sz = N/2, and is insensitive to total atom num- 
ber (Fig. [IJ . The atom-resonator coupling also gives rise 
to a differential light shift between the clock states, which 
we use to verify experimentally the coupling strength cal- 
culated from first principles from spectroscopically deter- 
mined resonator parameters. We measure a phase shift 
of 250(20) jirad per transmitted photon for a maximally 
coupled atom (on the resonator axis at an antinode of 
the probe standing wave), in excellent agreement with 
the calculated value 253(8) jirad [21 1. 

To account for the spatial variation in coupling be- 
tween standing-wave probe light and atoms, we define an 
effective atom number Nq — {{ri)'l/{if)e)Na, ~ 0.667Va, 
where the single-atom cooperativity 77, proportional to 
the local intensity of probe light, is averaged over the en- 
semble containing TVa atoms 2l|. The definition is cho- 
sen so that the projection noise variance of the effective 
atom number measured via the mode shift uj oc Na{ri)e 
satisfies the usual relation AA^q = Nq. This avoids carry- 
ing near-unity factors through the equations and allows 
direct comparison to a spatially uniform system of collec- 
tive cooperativity Ngijes, where ijes = (2/3)(r?^)e/(7?)e = 
0.47(1)770, taking into account the oscillator strength 2/3 
of the D2 line and the measured rms transverse cloud 
radius of 8.1(8) pm ^ w. The mode frequency shift per 
effective atom of population difference N between the 
clock states is dio/dN = 4.5(2) x 10"^ k [2lj. 

To quantify spin squeezing we need to measure AS^ 
and |(S)|. The latter can be obtained from the ob- 
served contrast C of Rabi oscillations as |(S)| — CSo, 
where the maximum spin Sq = No/2 is measured by 
optically pumping the atoms between the two hyper- 
fine states F = 1,2. For large Sq the cavity shift uj 



exceeds k {u! < 1.8k), which we take into account by 
correcting for the (accurately measured) Lorentzian line- 
shape of the resonator. To verify the atom numbers 2S'o 
thus obtained, we have also directly measured the cav- 
ity mode frequency shift w oc 5*2, finding agreement to 
within 2(4)% [2l|. AS'f is obtained from transmission 
measurements that always remain in the linear regime, 
with 2\ASz\duj/dN < 0.01k. 

The probe laser is frequency-stabilized to a far de- 
tuned, negligibly shifted mode [ilj. Each measurement 
of Sz employs two probe light pulses of duration T = 
50 jis = 158 ns separated by a 280 ps delaVjdur- 

ing which we apply a microwave tt pulse sequence [2l| to 
suppress inhomogeneous light shifts (spin echo sequence) . 
Each probe light pulse contains 10^ to 10^ photons which, 
after traversing the resonator, are detected with an over- 
all quantum efficiency Qe = 0.43(4). From the detected 
photon numbers in the two pulses, we deduce two cavity 
shifts uj± whose difference constitutes a single measure- 
ment M of S'z = {uj+ — uj-)/{4:duj/dN). In a typical ex- 
periment (Fig. 1(c)), after initializing the ensemble spin 
state by optical pumping into |1) (A) and applying a 7r/2 
microwave pulse to rotate the CSS into an equal superpo- 
sition of |1) and |2) (B), we perform two measurements 
Ml and M2 to induce and verify conditional spin squeez- 
ing. We quantify spin noise ASz by extracting variances 
from 100 repetitions of such a sequence. 

We determine the CSS projection noise level AS^^gg = 
Nq/A from the measured atom number A'o and verify it 



|6|, |15|, |17[ either by evaluating the variance Var(Mi) 
of the set of single measurements Mi; or by inserting 
between two measurements Mi and M2 a second CSS 
preparation, consisting of optical pumping into state |1) 
and a 7r/2 pulse, and evaluating Var(Mi — M2)/2. Fig. [5] 
shows the dependence of the corresponding quantities in 
atom number units, yi = 4Var(Mi) (open triangles) or 
2/2 = 2Var(Afi — M2) (open circles), on Nq. The con- 
tribution of CSS projection noise scales as AS^^gg oc 
Nq, while atom- number-dependent technical noise, e.g. 
due to microwave power fluctuations or any sensitiv- 
ity to atom number fluctuations, generically scales as 
ASztcch ^ ^0 • A quadratic fit yi,2 = ao + cliNq + a2NQ 
yields ai = 1-3(1) and 02 = 1(2) x 10~^ (not shown in 
Fig. [2]), but the data are also well fit by setting ai = 1, as 
required by independently measured cavity and atomic 
properties with no free parameters 2l|, and allowing 



a small technical noise contribution a2iVg < A'o with 
0-2 = 9(3) x 10~^ (solid curve). Slow drifts in microwave 
power of 0.4% over the set of measurements could account 
for the technical noise of yi , which vanishes if the data are 
analyzed by comparing only adjacent cycles of the exper- 
iment [2l|. Our ability to prepare an unentangled state 
close to a CSS — with Sz variance ASl^^^^ ~ 1.35*0/2 for 
our largest atom number — is not a prerequisite for spin 
squeezing but does provide independent confirmation of 
the CSS reference level for spin noise measurements. We 
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FIG. 2: Spin noise measurements (see text): yi — 4Var(Afi) 
for a single CSS preparation (open triangles), 3/2 ~ 2Var(Mi — 
M2) for two independently prepared CSSs (open circles), and 
4AS'f[„oas = 2Var(Mi — M2) for two measurements after a 
single CSS preparation (solid diamonds), all in units of atom 
number. Vertical error bars are statistical; horizontal error 
bars indicate standard deviation of measured atom numbers. 
The solid (dashed) line corresponds to ai = l,a2 = 9 x 10~^ 
(a2 = 0). Inset: Variance AS^ of Sz after rotating a squeezed 
state by an angle a about (S), with parameter-free theory 
curve. 



emphasize that, in quantifying spin squeezing below, we 
conservatively normalize to the CSS noise 4AS'^(^gg = A^o 
as obtained from our cavity parameters (dashed line) , not 
to the 30% larger slope of the unconstrained quadratic 
fit to yi^2- 

To prepare a state with (conditionally) reduced AS*^ 
(Fig. l(c)C), we simply measure Sz for a CSS on the 
X-axis with a photon number p « 5 x 10^ sufficiently 
large to resolve Sz beyond the CSS variance. Each such 
measurement Mi yields a value of Sz that is random but 
known, as verified by a readout measurement M2. We 
plot 2Var(Mi — M2) vs. atom number Nq in Fig. [5] (solid 
diamonds), finding it a factor of 2 above the photocurrent 
noise level, with very weak dependence on atom number, 
and well below the CSS level. 

In principle it is possible for the value of Sz at the end 
of the measurement to differ from the average value of Sz 
during the measurement. Besides the far-detuned lock- 
ing light whose effect on Sz is negligible, only spin-echo 
microwave composite tt pulses, whose fidelity was sepa- 
rately measured to be 98(1)%, and probe light are applied 
during Mi . The probe light can only change Sz through 
free-space scattering, which adds at most 3.1(3)% of CSS 
5 X 10^ 



projection noise at p = 5 x 10° [13|, |21|. Thus, while the 
added noise is negligible compared to the CSS level, it can 
explain part of the small remaining variance of Mi — M2. 

Provided Mi does not alter the state appreciably, and 
the measurements Mi , M2 are identical and uncorrelated 



variance of the state after measurement Mi can then be 

shown to be [AS^zUh = ^S^^^^AS^zmeJi^S^zprcp + 

AS'i^^) WM- When no new information is gained 
in measurement Mi (AS'^^^^^^ > AS'^^^.^^), the vari- 
ance is that of the state preparation process, AS'^p^.^p = 
Var(Afi) — ASzj^cas {close to, but above, the CSS value), 
while information gained reduces the variance, ultimately 
to the measurement variance AS'^j^^^^g of Mi. At A^o = 
3.3(2) X 10* and p = 6 x 10^, we observe a normalized 
spin noise ct^ = [AS'2]mi /AS'^^^gg = -9.1(8) dB (see Fig. 
Bl); a slight correction for the effect of photon scattering 
[2lj yields cr^ = -8.8(8) dB. 

The reduction of [AS'^]a/i below the CSS value ASzc^g 
is accompanied by a substantial increase in ASy because 
the differential light shift of the atomic levels, corre- 
sponding to a rotation of the Bloch vector about the z 
axis, depends on the intracavity intensity, which in turn 
depends on Sz- To observe the antisqueezing, we apply a 
microwave pulse after the squeezing measurement (at X 
in Fig. 1(c)) to rotate the spin state by a variable angle 
a about (S) before reading out Sz- The variance AS*^ of 
Sz in the rotated state, displayed in the inset to Fig. [21 is 
a sinusoid that is well described with no free parameters 
by our model of the ensemble-cavity interaction 21 1 . 

To verify spin squeezing, we also need to measure 
|(S)|, observable as the interference contrast C = |(S)|/S'o 
of Rabi oscillations induced between measurements Mi 
and M2- Fig. [3] shows C as a function of photon 
number p used in the state-preparation measurement 
at A^o = 4.0(1) X 10"^, and we have verified that the 
contrast C is independent of atom number [2l[. Both 
normalized spin noise and C can be fit by simple 
models (dashed and dotted curves) [21]. From these 
two measurements, we deduce the metrological squeez- 
ing parameter (solid triangles and solid curve). For 
p = 3 X 10^ we achieve Cm^ = C^/icr^Cin) = 3.0(8) dB of 
metrological gain (and an inverse entanglement parame- 
ter Cc^ = C/cr2 = 4.2(8) dB, not shown). The finite ini- 
tial contrast Cjn = Si^/ Sq — 0.7 in the ensemble without 
squeezing is due to the resonator locking light, and can 
be improved by detuning this light further from atomic 
resonance. The probe-induced contrast reduction proba- 
bly arises from differential light shifts between the clock 
states that are imperfectly canceled by the spin echo tech- 
nique because of atomic motion. In the absence of any 
technical noise, a fundamental limit to the spin squeez- 
ing, associated with photon scattering into free space, 
would be Cm^ < ^y {3/2)Nor]eB ~ 18 dB in our system 



13|,ll4|,|21| 



21|, AS^^^^ = Var(Mi - M2)/2 represents the uncer- 



tainty of any single such measurement. The conditional 



with cooperativity iVoryoff ~ 3100 

For the data presented above, the readout quantifying 
the entanglement was completed 500 ps after prepara- 
tion of the squeezed state. We have further verified that 
the squeezing remains after a Ramsey clock sequence, in 
which two tt/2 pulses about the cc-axis, separated by a 
short (70 ps) precession time, are inserted at X in Fig.[TJ;. 
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FIG. 3: Measured data with fits of simple models ^2l| for nor- 
malized spin noise cr^ = [A5f]A/i/A5f qss (open diamonds), 
contrast C (open squares), and metrological squeezing param- 
eter (solid triangles). 



Such a clock can achieve precision below the SQL because 
the first of these 7r/2 rotations initiates it with a phase 
that is known, from the squeezing measurement, to bet- 
ter than the CSS uncertainty. 

The phase coherence time of the unsqueezed CSS in 
our current trap is 10(2) ms. Both microwave and op- 
tical clocks with 1 s coherence times have alr eady 
been demonstrated with trapped atoms 
Whether and to what degree the squeezing technique 
demonstrated here could benefit such clocks and other 
precision experiments Q will depend on the clock char- 
acteristics, noise sources [l6| . and lifetime of the squeezed 
state. These questions, as well as possible systematic ef- 
fects, need to be investigated in the future. 

The group of E. Polzik independently and simultane- 
ously achieved results similar to ours in a Mach-Zehnder 
interferometer [2^. We have recently demonstrated a 
new squeezing method using cavity feedback [30]. 
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I. OPTICAL RESONATOR AND DIPOLE TRAP 

The parameters of the near-confocal Fabry-Perot res- 
onator at the wavelengths of trap hght (851 nm) and 
probe hght (780 nm) are summarized in table lAll Both 
trap laser and probe laser are locked to the resonator via 
Pound-Drever-Hall locks of 1 MHz loop bandwidth. An- 
other feedback loop stabilizes the resonator frequency to 
an atomic transition in '^^Rb. 

The atoms are loaded into the standing-wave opti- 
cal trap from a microchip-based magnetic trap described 
elsewhere pj. After polarization gradient cooling in the 
linearly polarized optical trap, we apply a 5.6 G magnetic 
field along the resonator axis and a circular polarization 
fraction of 0.5(1) to the trap light. This combination 
yields a first-order cancellation of the vector and scalar 
light shifts, minimizing inhomogeneous broadening of the 
clock transition. Table IA2I summarizes the characteris- 
tics of the atomic cloud in the optical dipole trap. 



II. DETECTION SETUP 

We probe the atoms' index of refraction on the D2 
transition with linear polarization through the optical 
cavity. The probe laser carrier lies halfway between two 
TEMqo modes of the resonator. A broadband electro- 
optic modulator (model PM-0K5-10-PFA-PFA-780-UL 
from EOSPACE) is used to generate sidebands for lock- 
ing and probing (see Fig. IA1[) . A lock sideband at 



Parameter 


A = 780 nm 


A = 851 nm 


Mirror separation L 


26.62(1) mm 


Mirror curvature radius R 


25.04(2) mm 


Free spectral range (jJfskJ (Ztt) 


5632.0(2) MHz 


Transverse mode spacing u)t/(2Ti) 


226.3(3) MHz 


Linewidth k.x/{2-k) 


1.01(3) MHz 


135(2) kHz 


Finesse Tx 


5.6(2) X 10^ 


4.2(1) X 10^ 


Mode waist wx 


56.9(4)iim 


59.5(5)iim 


Antinode cooperativity rjo^x 


0.203(7) 


1.65(4) 



TABLE Al: Resonator parameters. The mode waists are 
calculated at the position of the atoms. Outside this table, 
all resonator values refer to the probe wavelength A = 780 nm. 



113 MHz, resonant with a TEMio mode, produces the 
Pound-Drever-Hall error signal. 

The probe sideband at (5wfsr + k)/2 « 27r x 
14080 MHz lies on the slope of a TEMqo resonance with 
a detuning of ~\-2tt x 3.57(1) GHz relative to the atomic 
F = 2 ^ F' = 3 transition. The far off-resonant sym- 
metric (compensation) sideband at — (SwpsR + '*)/2 lies 
on the opposite slope of another TEMqo mode, such that 
the total transmission in the two modes is (ideally) sensi- 
tive only to atom-induced shifts of the cavity resonance, 
but not to frequency jitter of the laser relative to the 
cavity. 

The transmitted power in the TEMqq mode is cou- 
pled into a single-mode fiber to filter out the lock side- 
band and subsequently detected with overall quantum 
efficiency Qe — 0.43(4) on a Si avalanche photodiode 
(Hamamatsu model S3884). At a typical power of 2 nW 
for our T = 50 ps long probe pulses, the total photodetec- 
tion noise (including excess noise of the avalanche pho- 
todiode operated at gain M = 13 0) is a factor of 1.9 in 
variance above the photocurrent shot noise. 



III. STATE PREPARATION 

Each cycle of the experiment includes three CSS prepa- 
rations with the same loaded atoms. The first CSS prepa- 
ration precedes the measurements Mi and M2 used to 
prepare and read out a squeezed state. The measurement 
All is also used to quantify the unconditional variance of 
Sz via j/i = 4Var(Afi). The two subsequent CSS prepa- 
rations precede the measurements Mi and M2 used for 
independent verification of the state preparation noise 
via 2/2 = 2Var(Afi - A'h)- 

First CSS preparation: The atoms are optically 
pumped into |1,0) using (T-|_/(T_ -polarized light on the 
F = 1 ^ F' = Q transition while repumping ox\ F — 
2 — > F' = 2. To improve the state purity, a (composite 
SCROFULOUS [3]) microwave tt pulse is subsequently 
applied on the |1,0) — > |2,0) transition, all F = 1 states 
are emptied using resonant light on the F = 1 — >• F' = 1 
transition, a second tt pulse returns atoms from |2, 0) to 
1,0), and all atoms remaining in F = 2 (^ 12% of the 
initial atom number) are expelled from the trap using 
resonant light on the F = 2 — F' = 3 transition. After 
this procedure, more than 99% of the remaining atoms 
are in the state |1, 0). A 7r/2 pulse prepares the atoms in 
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Optical dipole trap 


Atomic cloud 


Axial frequency 




550 kHz 


Length 


I 


1 mm (~ 2000 wells) 


Radial frequency 


CJr/(27r) 


1.8 kHz 


RMS radius 


cr,. 


8.1(8) iim 


Trap Depth 


Uo/h 


24(1) MHz 


Radial temperature 


kBT,/h 


1.3(2) MHz 



TABLE A2: Characteristics of standing-wave dipole trap and atom cloud. The trap depth and trap frequencies are determined 
from the intracavity power and mode geometry. The radial temperature is measured by suddenly releasing the atoms and 
observing their ballistic radial expansion as a decrease in coupling to the resonator. 
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FIG. Al: Laser stabilization and detection scheme, indicating frequencies of carrier and lock sidebands (red) and probe and 
compensation sidebands (blue) relative to cavity resonances and atomic transitions. Not to scale. 



a CSS with {Sz) = 0, on which the squeezing and readout 
measurements are performed. 

Second and third CSS preparations: We then pro- 
ceed, using the same loaded atoms, to prepare a CSS in 
the xy-pleme two more times to confirm the CSS projec- 
tion noise. In order to compare two identically-prepared 
CSSs with the same total atom number, we forego the 
state purification procedure described above, since it 
leads to a ^ 12% loss of atoms. We thereby allow our 
imperfect optical pumping to leave 12(2) % of the atoms 
in |1,±1). These residual atoms do not contribute to 
our measurement of Sz because our spin echo technique 
(see Sec. |vi| cancels any contribution from atoms not ad- 
dressed by microwaves resonant with the |1,0) — >■ |2, 0) 
transition. Therefore, in Fig. 2, the atom number A^o 
for the data derived from this pair of preparations (open 
circles) includes only atoms in |1, 0) and |2, 0) and is sys- 
tematically 12% lower than the atom number used to 
measure Var(Mi) and ^S^^^^^ (open triangles and sohd 
diamonds). 



State Preparation Noise 

Figure 2 indicates the presence of technical noise in 
the state preparation. The technical noise evident in 
Var(Mi) is probably due to slow drifts in microwave 
power. An alternative analysis, in which we compare 
each measurement Mi with the value Mf^""^ in the pre- 
ceding experiment cycle, yields a fit 2Var(Mi — Mf"^°'^) — 
2650(400) + 0.95(23)Aro + 1(9) x IQ-^N^, i.e. a result 
consistent with no contribution from the quadratic term. 
We also show in Fig. 2 the variance y2 — 2Var(M2 — Mi), 
which is immune to slow drifts in microwave power. How- 
ever, after completing this work, we discovered that the 
state preparation preceding the measurement M2 was 



compromised by an effect of leakage light during that 
preparation, to which we attribute the small technical 
noise observed in 2/2- 

IV. ATOM-LIGHT INTERACTION IN AN 
OPTICAL RESONATOR 

We summarize the theory of the interaction of a two- 
level atom with an optical resonator mode at large detun- 
ing (5 ^ r, relative to the excited-state linewidth F, from 
the atomic transition. The extension to our real system 
of many atoms with nontrivial level structure follows in 
Sec. m 

A. Atom-Resonator Coupling, Cooperativity, and 
Optical Depth 

The atom-resonator coupling g{r) = \deg ■ E(r)|/fi, for 
an atom at position r — (p, z) in the Gaussian mode is 
given by 

where dgg is the dipole matrix element between the two 
states |f;) and |e), ujeg is the energy of the transition, w{z) 
is the mode waist at the position of the atom, and L is 
the resonator length. {2g is the vacuum Rabi frequency.) 
The coupling g{r) is related to the atomic excited-state 
linewidth F — uj^gd^g / {Sireahc^) and resonator linewidth 
K by the single-atom cooperativity ri(r), the ratio of the 
scattering rate into the resonator mode to the free-space 
scattering rate Q: 

.(r) = ^ = ^e-^V»%i,.(,,), (2) 
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where T = ttc/ILk) is the finesse of the resonator and 
k = uJeg/c is the probe wavenumber. 

The cooperativity is closely related to the resonant op- 
tical depth, which for a single atom with scattering cross 
section CTsc in a uniform beam of area A in free space is 
given by asc/A. A light pulse resonant with the cavity 
passes through the atomic sample on average 2F / tt times. 
For an atom at an antinode of the standing-wave mode, 
the resonator then enhances the resonant optical depth 
by a factor of AT/tt relative to its value 12/(fc^w^) on 
the axis of a running-wave Gaussian beam of waist w in 
free space, so that 2r] represents the resonator-enhanced 
single-atom optical depth. 



Resonator Mode Shift and Back-Action Phase 
Shift 



iVo ~ l^a, is chosen such that the projection noise vari- 
ance satisfies the usual condition for a uniform sample. 



No 



(4) 



The mode shift due to an ensemble in state \F) is 
^{P) — NQrics^K/i^Sp), where ■q^e is related to the co- 
operativity 770 = 24J^/(7rfc^?i;^) of a maximally coupled 
atom by 

^ / / . 2 , \ 2 I g 2 ^ 0.47(1) (5) 

for our cloud of radius Ur w. 



Measurement of Sz 



We now consider a resonator containing n photons 
and a single atom in state \g). The shift oji — g^/S = 
•qTn/^iS) of the resonator mode due to the interaction 
with the atom is accompanied by an AC Stark shift nwi 
of the atomic level \g) due to the light; the symmetry 
between these two effects is readily understood in the 
dressed-atom picture Since photons are transmitted 
through the resonator (leaving the resonator in the for- 
ward direction) at a rate nK/2, the phase shift of the 
atomic state \g) per transmitted photon is 2uji/k. 



V. POPULATION MEASUREMENT 
A. Mode Shift and Effective Atom Number 

In Sec. Hvl we expressed the shift loi = g^ /5 of a res- 
onator mode coupled to a two-level atom in terms of 
the cooperativity rj. In any real atom, at finite detun- 
ing, the coupling g is polarization-dependent and must 
be summed over various excited states. In terms of the 
cooperativity 77 of a two-level atom (i.e. the cooperativity 
on a cycling transition), given by the right-hand side of 
Eq. [21 a single *^Rb atom at position r occupying state 
\F) shifts the mode frequency for linearly polarized light 
on the D2 transition by an amount 



The probe sideband is tuned to the frequency between 
the 15251/2, - 1) ^ 152^3/2) and \^^S^/2.F = 2) ^ 
15^^3/2) transitions at which the atom- induced differen- 
tial frequency shift between the probe and compensation 
sidebands, 



{N2 - Nl)r^,sT^^ 
4(5' 



(6) 



(r)=/^(r)^ 



(3) 



is proportional to the effective-atom population differ- 
ence N = N2 — Ni between the hyperfine states F — 1,2, 
but independent of the total atom number A^i -I- 7V2. 
Here S' = 2tt x 3200(10) MHz; see Table [A3] for de- 
tails. For our cloud geometry and probe polarization, 
the differential mode shift per effective atom is dcu / dN = 
2n X 45(1) Hz/atom = 4.5(2) x 10"^ K/atom. The mode 
shifts due to projection noise on N are much smaller 
than k/2, leaving the resonator transmission in the lin- 
ear regime, such that the change in transmitted power 
is directly proportional to Sz = N/2. Nevertheless, in 
all measurements of Sz we take into account the full 
Lorentzian lineshape of the resonator transmission, al- 
lowing the same procedure to be used for measuring pro- 
jection noise, contrast, and total atom number. The cav- 
ity linewidth k/{2tt) — 1.01(3) MHz at the probe wave- 
length is accurately measured by tuning the probe side- 
band over the TEMqo resonance of the empty cavity and 
measuring probe transmission. 



Here, T is the excited-state linewidth; Sp is an effective 
detuning from the 15^51/2, i^) \^'^Ps/2,F') transitions 
averaged over excited hyperfine states F'; and / = | 
is the oscillator strength of the D2 line. For our ex- 
tended sample of A^a atoms we define the effective coop- 
erativity 77eff = f{'r]^)e/{r])e and the effective atom num- 
ber A^o = Na{r])l/ {r]'^)e- Here, ()e denotes an average 
over the atomic ensemble. This definition, which yields 



Measurement of N( 







At the end of each experiment cycle, we measure the 
effective atom number A^o- We determine the atom num- 
ber by pumping all atoms first into F = 2, then into 
F = I, and in each case measuring the resonator trans- 
mission. Although the resonator mode shifts w^^^ are 
linear in Nq, they are on the order of k for our typical 



4 



Resonance 


<523/(27r GHz) 


Cavity Shift per Effective Atom/(K x 10"®) 


\F = l,mF = 0) 


|F = 2,mF = 0) 


Probe 
Lock 

Compensation 


3.57(1) 
-10.40(1) 
-24.59(1) 


-49(2) 
-0.6(1) 
-4.6(2) 


-1-39(1) 
-1.0(2) 
-5.8(2) 



TABLE A3: Frequency shifts per effective atom in either of the clock states for relevant resonator modes. ^23 is the detuning of 
the given mode from the \5'^Si/2,F = 2) — > P3/2, F' = 3) atomic transition. The calculated shifts include the excited-state 
hyperfine structure, as well as the spatial overlap of the cloud with the mode. 



atom numbers, so that the transmitted power is a non- 
hnear function of w^^^ which we must invert to obtain 
the atom number. To verify our determination of the 
resonator shift from the non-hnear transmission signal, 
we additionally measure the average mode shift over sev- 
eral cycles of the experiment by finding the probe side- 
band frequencies ujp that maximize the probe trans- 
mission when all atoms are pumped into hyperfine state 
F. By thus directly measuring the mode frequency shift, 
we obtain a linear measure of effective atom number 
ATjf = (wf ^ - 4^V2/(27r X 45 Hz). In Figure |A2l 
we plot Nq against the average atom number Nq ex- 
tracted from the nonlinear transmission signal. The fit 
= 0.98(4)iVo with reduced = 0.3 indicates that 
the two measurements are in good agreement. 

At the large atom number A^o = 3.3(2) x 10^ where 
we calculate squeezing parameters, the dominant uncer- 
tainty in iVo arises from sensitivity to the initial place- 
ment of the probe and compensation sidebands. In ex- 
tracting A^o from the transmission, we assume that we 
have correctly placed the probe and compensation side- 
bands at ±k/2 detuning from cavity resonance when 
Sz = 0. We monitor the placement of the sidebands 
in each cycle of the experiment by shifting the frequency 
of the laser relative to the resonator by +k and —k from 
the usual configuration and measuring the transmission 
in each case. Hence, we are confident that any systematic 
error in the placement of these sidebands is less than the 
shot-to-shot fluctuations. We therefore always estimate 
the uncertainty in A^o by the standard deviation of the 
calculated Nq values. 



D. Experimental Verification of Atom- Resonator 
Interaction 

The calculated mode shift per atom, Eqs. [3] and 
[6l is used to convert measured transmission into atom 
number. We verify it by measuring the complementary 
atomic phase shift cj) = 2(wp' -~lo[^'^)/k induced between 
states |1) and [2) by a single probe photon transmitted 
through the resonator. We determine the phase shift (po 
of a maximally-coupled atom by means of a Ramsey mea- 
surement 6] , applying an optical probe pulse of variable 



35x1 oV 





30- 






CD 

E 


25- 


CD 




13 

cn 


20- 


c« 




CD 

E 


15- 






m 

CD 


10- 


c 






5- 


_i 








Ot 



5 10 15 20 25 30x10 
Nq (non-linear measurement) 

FIG. A2: Comparison of two methods of determining effective 
atom number, as described in the text. The vertical error bars 
indicate the uncertainty in positioning the probe sideband on 
cavity resonance. The horizontal error bars arise because the 
atom number No is the average over a different set of loading 
cycles from those in which Nq is measured. 
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FIG. A3: Measurement of the atom-photon interaction. A 
probe pulse of varying photon number p is inserted into a 
Ramsey sequence, resulting in a differential light shift between 
atomic states. The light-induced phase shift and decoherence 
(open circles) can be suppressed by a spin echo technique 
where a microwave n pulse is inserted between two probe 
pulses (solid squares). 



5 



duration between two microwave 7r/2 pulses. The pop- 
ulation difference 2Sz, measured via the resonator shift 
(Fig. IA3p . is an oscillatory function of the transmitted 
probe photon number p. The oscillation is damped due 
to inhomogeneous light shifts. For an ensemble of atoms 
on the resonator axis evenly distributed with respect to 
the probe standing wave, a spin state prepared along the 
X-axis of the Bloch sphere acquires, after the interaction, 

Jq cos (p0o sin^(fcz)) sin^(fcz) 
sin (kz) dz 
=Jo (u) cos (u) — Ji (u) sin {u) , (7) 

where the J„ are Bessel functions of the first kind 
and u — p<j)^/2. From a fit of this form we extract 
4>o = 230(20) jirad. (The phase shifts due to lock 
and compensation light are negligible.) A fit to a full 
numerical model including the radial cloud size yields 
0moas _ 250(20) jirad, in excellent agreement with the 
value (/ig^''^ = 253(8) prad calculated from cavity param- 
eters. 

VI. SPIN ECHO SEQUENCE 

We use a spin echo technique to reduce the probe- 
induced inhomogeneous broadening (Fig. IA3I) . All 
probe light is applied in two 50 ps long pulses sepa- 
rated by a composite tt pulse, consisting of a sequence 
i?^/3(7r)i?_^/3(7r)i?^/3(7r) of three simple microwave tt 
pulses, where the subscripts indicate phases chosen to 
compensate variations in pulse area Q . The spin echo is 
optimized at a probe pulse separation of 330(20) ps, cor- 
responding to a half-period of the radial trap oscillation. 

VII. DATA ANALYSIS 
A. Noise Model 

The spin measurement variance ^S^^^,^^ — Var(Mi — 
M2)/2 as a function of probe photon number p used 
in the measurement (Fig. IA4[) is well described by as- 
suming the following independent noise contributions to 
4AS'2jj^^^^ = 6_2P"^ + b^ip-^ + 6oP° + bip^: electronic 
noise of the detector scahng as ^S-^^^^ cx pho- 
ton shot noise and avalanche excess noise scaling as 
^•^zmcas ^ P~^', technical noise that is independent of 
photon number, including the effects of imperfect mi- 
crowave rotations in the spin echo procedure; and noise 
due to photon (Raman) scattering, AS'^^^^^ oc p. We 
quantify each of the known noise contributions: 

• Photon shot noise and avalanche excess 
noise: We calculate the photocurrent noise due 
to the probe and compensation light, taking into 
account both the shot noise of the light detected 




10' 10' 
Probe photon number 

FIG. A4: Comparison of observed variances with noise model. 
Shown as a function of probe photon number p are the 
measurement variance in atom number units, iAS'^^^^^ = 
2Var(Afi — M2), with (solid blue diamonds) and without 
(open green squares) atoms; and the variance 4Var(Mi) with 
atoms, which includes projection noise. In addition, vari- 
ances 2Var(Mi+ — M2+) (solid gray inverted triangles) and 
2Var(Mi_ —M2-) (solid purple triangles) using a single probe 
pulse from each measurement (as shown in Fig. IA5|) are plot- 
ted as a function of single-pulse probe photon number p/2. 
Curves correspond to the noise model described in Sec. 1X1 

with a quantum efficiency Qe = 0.43(4) and the 
excess noise factor /apd = 1.9(4) of the avalanche 
photodiode operated at a gain of 13 0], obtaining 

b-i = 2(/APD/Qe)(rfiV/rf(2w/K))2 = 1.1(3) X 10-9. 

• Electronic noise: From a fit to the noise mea- 
sured in the absence of atoms (open green squares 
in Fig. lA4p in which we constrain the coefficient 5_i 
to the value calculated above, we obtain an elec- 
tronic noise contribution h-2 = 6(1) x lO^'^, most 
of which is attributable to the Johnson noise of the 
transimpedance gain resistor in the photodctection 
circuit. 

• Microwave infidelity: The tt pulse used in the 
spin echo produces at least 98(1)% inversion. We 
model the imperfect tt pulse as a perfect one com- 
bined with an incoherent process that flips on av- 
erage ^ — 2(1)% of the spins, yielding 60, ^ = 
^jlNq. We treat the errors as incoherent because 
the atomic phase is inhomogeneously broadened 
by {4>o/2)p = 1.3 X IQ^^p radians when the mi- 
crowaves are applied. At the optimum photon num- 
ber p = 3 X 10^ for squeezing, the Ramsey contrast 
remaining after a single probe pulse is only 10(3)% 
(see Fig. IA3I) . We briefly address possible coherent 
noise processes below (Sec. placing an upper 
bound on the effect of such processes. 

• Raman scattering: In our system, the prob- 
ability of a Raman scattering event is PRam = 
5.6 X 10~^ per probe photon transmitted through 
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the resonator. This value, calculated including the 
full excited-state and ground-state hyperfine struc- 
tures, includes probabilities PRam = PAF+PAm,F + 
PafAwlf corresponding to three types of scatter- 
ing events: those which change F but not mp, 
those which change mp but not and those which 
change both F and m^?, respectively. To first or- 
der in these probabilities, the total contribution of 
Raman scattering to the measurement variance is 

bl = (4/3PAF + l/2FAm^ + l/iPAFAn^M = 

4.7 X 10~^-/Vo per probe photon. Section 1X1 outlines 
the derivation of this expression for hi . 



Figure IA4I shows a fit of the above model (solid blue 
curve) to the observed measurement variance in atom 



number units, ^^Sj. 



at TVo = 3.3(2) x lO" and vari- 



able photon number p (solid blue diamonds) . We fix the 
noise contributions enumerated above and leave free a 
term 6o,tec?i to account for technical noise that is indepen- 
dent of probe photon number but not due to microwave 
infidelity. The value bo,tech = 1400(400) = 0.04(l)Afo ob- 
tained from the fit may be due to frequency jitter of the 
laser relative to the cavity that is incompletely canceled 
by the compensation sideband; or to technical noise in 
the probe light level, e.g. from fluctuations in the cou- 
pling of the light to the cavity. The data show very good 
agreement with our noise model with this single free pa- 
rameter bo^tech, all other parameters being independently 
measured or calculated as described above. 

The distinction between the two contributions to &o = 
^o,tcch + ^o,Ai is verified by also plotting measurement vari- 
ances (in atom number units) 2Var(Mi_ — M2-) and 
2Var(Mi+ — M2+) obtained by comparing either the first 
and last or the second and third of the four probe pulses 
Mi± constituting the two measurements Mi, i £ {1,2}; 
see Fig. IA5I Whereas two tt pulses are applied between 
Ml- and M2^, no microwaves are applied between Mi+ 
and A/2+- Thus, fits (not pictured) to 2Var(Mi_ - M2_) 
and 2Var(i\fi+ — M2+) reveal the microwave infidelity 
;^fit = (&- - 6+)/(4iVo) = 0.03(1), consistent with the 
independently determined value = 0.02(1). (The value 
/j,fit also includes a small contribution < 0.01 due to pho- 
tons scattered by the lock light.) 

In addition, we plot 4Var(Mi) (open red circles in 
Fig. IA4p and a curve given by the following expres- 
sion: 4Var(Mi) = + b_i/p + &o,tech/prep + (1 - 

A* - [2PAF/3 + PAm^/2 + 2PAf Am^/3]p)iVo (sec Sec. El- 
Here, we fit the term 6o,tcch/prcp = 0.14(7)iVo to allow for 
technical noise in the state preparation but constrain all 
other parameters to the values given above. 

The red curve in Fig. 3 is derived by combining the fits 
to 4AS'^j^gg^g and 4Var(Mi) described above with a fit to 
the contrast, described in Sec. |d1 in accordance with Eq. 
[T5] below. The dashed blue curve in Fig. 3 corresponds 
to the same expression without the contrast factors. 
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FIG. A5: Measurements Mi and A/2 of Sz for spin squeezing 
and readout. Each measurement is obtained from the probe 
light transmitted (solid blue line) in two pulses, before and 
after applying microwaves (dashed red line) to perform a tt 
rotation. Specifically, each measurement can be expressed as 
the average Mi — {Mi+ + Mi-)/2 of two single-pulse mea- 
surements Mi± = ±u]i±/\2dui / dN\, where u]i± are the atom- 
induced shifts of the cavity resonance deduced from the trans- 
mitted probe light. With each vf rotation, we switch the sign 
convention for Sz = ±(-/V2 — Ni)/2 to compensate for the 
population exchange between states |1) and |2). 



Bound on Coherent Microwave Errors 



We have noted above that any noise process associated 
with coherent microwave pulse errors in the spin echo 
is highly suppressed by inhomogeneous broadening of 
the atomic phase before the spin-echo- induced rephasing. 
For example, the drifts of 0.4% in microwave power (0.2% 
in Rabi frequency) which we infer from our state prepa- 
ration noise AS'^p^.^^p would lead to rotation errors sub- 
stantially less than i50niax = 2 x IO^'^tt in our composite 
TT pulse, which is designed to compensate for microwave 
power errors. At the optimum squeezing photon number 
p = 3 x 10^, where the interference contrast is Cse ~ 10% 
when the spin-echo microwaves are applied (see Fig. IA3p . 
the resulting normalized variance in would be at most 
A52^„h/A^fcss < S(l>l,ANo = 4 x IQ-^No « 0.01 
at No 

FA c2i_ /At, 

photon number. 



33000. This is much smaller than the normal- 
ized spin noise [AS'^]mi / ^S^^^g = 0.20(3) at the same 



Note further that, while the calculation of condi- 
tional spin noise from the spin measurement variance 
A'^zmcas = Var(Mi — M2)/2 obscures measurement er- 
rors that are perfectly correlated, it overestimates by a 
factor of two the effect of measurement errors that are 
perfectly anticorrelated. Since correlated microwave er- 
rors can yield either correlated or anticorrelated measure- 
ment errors depending upon the atomic phase, which in 
turn depends sensitively on the number of probe photons 
applied before the tt pulse in each measurement, we ex- 
pect a further suppression of any correlated noise below 
the bound given here. 
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Derivation of Spin-Flip Terms 

This section outlines our derivation of the contribu- 
tions of incoherent spin flips to the variances 4:AS'^^^^^ 
and 4Var(Mi) in our noise model above. Both 
= Var(Mi+ + Mi_ - M2+ - M2_)/2 and 
4Var(Mi) = Var(Afi_|_ + are composed of covari- 

ances Gov (Mia, Mj^) of single-pulse measurement out- 
comes, where i,j G {1,2} and a, (3 G {+, — } (see Fig. 
IA5|) . To evaluate these covariances, we first express each 
single-pulse measurement M^q, as 



,{t)dt 



2\duj/dN\T 



(8) 



where cj^c (t) is the atom- induced cavity shift at a time t 
from the beginning of pulse ia; we here neglect errors in 
determining the atom-induced cavity shift due to photon 
shot noise and technical noise, as these are uncorrelated 
with the atomic state and can be treated separately. In 
the absence of Raman scattering, each ujia is constant in 
time; and if the spin-echo tt pulses have perfect fidelity, 
then uji^ = —uJi^. Raman scattering and microwave in- 
fidelity cause deviations from this ideal behavior, so that 



Cov{uia{t),ujji3{t')) « a/3 



duj 



dN 



7Vo[l-2r,„j^(i,i')] 



(9) 

Here, for a = /3, we define riaji3{t,t') as the proba- 
bility for an atom to be in a different hyperfine state 
F at time t' of pulse j/3 than at time t of pulse ia; 
for a — — /?, meaning that there is exactly one mi- 
crowave pulse between t and t' , riaj/3(t,t') is the prob- 
ability for the hyperfine state to be the same at both 
times. We have made the approximations that the ini- 
tial state preparation is projection-noise limited; and 
that riajp(t,t') <C 1, meaning that all the incoherent- 
spin-flip probabilities — namely, the Raman scattering 
probabilities pPaf, pPAmp, and pPAFAmp and the mi- 
crowave infidelity p, — are small. Each expression for 
riaj/3{t,t') is then, to lowest order, a linear combina- 
tion of these probabilities with coefficients that depend 
on how many tt pulses are applied between pulse ia 



and pulse j/S. For example, ri-^i^(t,t') « ^i^pP^F 

2rp-pPAFAmp 



\PAr. 



/i. This expression and simi- 



lar expressions for the other riaji3(t,t') are used to eval- 
uate the covariances Cov {Mia, Mj 13) using Eqs. Eland 
m We thus obtain the terms (4/3Paf + l/2PAnij, + 
l/iPAFAmp + n)No in our model for A/\Sl^^^^ and the 
terms [l ~ pi - [2Paf/3 + Pa^J^ + 2Fafa™f/3]p)A^o 
in our model for 4Var(Afi). 

These results can be understood qualitatively as fol- 
lows. Incoherent spin fiips increase the spin measure- 
ment variance l^S^^^^ = Var(Mi — Ah)/^ by reduc- 
ing the correlation between the squeezing and readout 
measurements Mi and M2. However, they diminish the 



projection noise observed in the variance Var(Mi) of a 
single measurement, since any scrambling of spins allows 
the measurement to average over the different ensemble 
spin states at different times. 



B. Derivation of Metrological Squeezing Parameter 

In interpreting ASl^^^^ as a measurement uncertainty, 
and in deriving the quantum uncertainty [AS'^J^v/j of the 
state prepared by the squeezing measurement, we make 
assumptions of uncorrelated noise in the two measure- 
ments Ml and M2. These assumptions are justified by 
the noise model in Sec. El The dominant noise contri- 
butions at the optimum squeezing point — photodetector 
noise and technical noise attributable to frequency jitter 
between laser and cavity — are all uncorrelated between 
the measurements. (Any frequency shaking that is slow 
enough to be common to both Mi and M2 is also com- 
mon to the cavity shifts ujj^ and a;_ before and after the 
spin echo and thus does not affect our measurement of 
Sz oc w+ — UJ-.) The remaining noise, due to changes in 
the atomic state via Raman scattering or microwave tt 
pulse infidelity, is correlated with the atomic state and 
hence affects the correlation between the two measure- 
ment outcomes. We here generalize the derivation of the 
conditional spin noise and metrological squeezing param- 
eter to encompass the small effects of such spin-flip noise, 
thereby fully accounting for all processes in our noise 
model. 

The goal of our analysis is to evaluate (^m = 
2[AS'^]Mi<S'in/(KS)P) using results of our measurements 
Ml and M2. Since the measurements can change the 
atomic state, we define Szf as the value of Sz at the end 
of the first measurement Mi . The conditional quantum 
uncertainty [AS'^Jmi ol t-^ic state prepared by the first 
measurement is found by minimizing Yax{Szf — wMi) 
with respect to the weight w given to the measurement 
information 



[AS! 



I Ml = mm, 



;{Var(5. 



^Mi)}. 



(10) 



The minimum occurs at w — Cov(iS'z/, Mi)/Var(Mi) and 
is given by 



[As^Uh = y^Szf) 1 1 



Cov(5,j,Mi)^ 
Var(5^/)Var(Mi) 



(11) 



As described in Sec. El the noise in our system is 
well described by noise that is uncorrelated with Szf and 
by incoherent spin flips. We begin with the simplify- 
ing assumption that all Raman scattering events transfer 
atoms between the two clock states. In this model, for 
identical measurements AIi and M2, 



Cov(5,/,Mi)2 = Cov(%,Mi)Cov(%,M2) 
= Cov(Mi,M2)Var(5,/), 



(12) 
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since the only correlations between the two measure- 
ments are due to their hnear dependence on the com- 
mon value Szf of Sz at the end of the first measure- 
ment and the beginning of the second measurement. 
Here, Yax{Szf) represents the unconditional variance of 
an ideal, noiseless readout following the squeezing mea- 
surement. 

We would like to express Var(S'z/) in terms of mea- 
sured quantities. We define ^Sl^^^ = Var(Mi - M2)/2 



and A5'2p^^p 



Var(Mi) - A5i,„_ = Cov(Mi,M2), 
where the latter equality holds because Var(Mi) = 
Var(Af2). In the limit where the measurement does not 
change the atomic state and all measurement noise is 
uncorrelated, A5^pj.j,p represents the unconditional vari- 
ance of Sz, including the CSS projection noise and any 
technical noise in the initial state preparation. To allow 
for a measurement that induces spin flips, we use Eq. [T^] 
to relate AS'^p^.^p to Y'Ar{Szf), obtaining 



^^fprep ^ / Cov(Af„5,/) y ^ 



(13) 



To lowest order in the spin-flip probabilities associated 
with Raman scattering and microwave pulse infidelity, 
tp « pPaf + M- The measured unconditional variance 
ASzp^cp lower than the true variance Var(S'2/) by the 
factor (1 — ep)^ because those spins that are flipped part- 
way through the first measurement contribute less to 
the measurement (of the time-averaged spin) than to the 
value Szf at the end of the measurement. Hence, 



A 92 AC2 

2 mcas 2 prep 



[l-e^nAS^zprep + ^Sl 



(14) 



The factor (1 — ep)^ has a 0.3 dB effect, smaller than the 
statistical uncertainty, on the minimum normalized spin 
noise = [AS^]mi/ AS^zcss^ at p = 6.4(6) x 10^ We do 
not include this factor in Fig. 3, but we correct for it in 
reporting the minimum normalized spin noise , as well 
as the entanglement parameter Co at p = 3 x 10^. We 
show below that the metrological squeezing parameter 
Cm is independent of Cp. 

Note that if we define the measurement strength 



AS* 



2: prep 



/ASz 



g, then for a projection- noise- 



limited state preparation (AS*^ 



^'S'^css) and a 



measurement that does not change the atomic state {cp 
0) Eq. [13] reduces to the_expression cr^ = 1/(1-1- n'iicas) 
given by Appel et al. |9]. Whereas Appel et al. de- 
fine the measurement strength for a photon-shot-noise- 
limited measurement, we have here extended the defini- 
tion to include more general noise. 

We wish to determine the metrological squeezing pa- 
rameter Cm = 2[AS'2]^^5'i„/(|(s)|2) where 1(S)| is the 
length of the mean spin vector at the end of the squeez- 
ing measurement. The Rabi oscillation curve from 
which we obtain the contrast (see Sec. |D]) is mea- 



sured using our standard readout, identical to the squeez- 
ing measurement. The contrast C{p) = |(S) 1/5*0 that 
would be observed in an ideal readout of Szf is related 
to our measured contrast Cmcas(p) by C^casiv) /C{p) = 
Cov(M2, S'^/)/Var(5^/) = 1 - Cp. Hence, 



Cin \ASl 



Ml 



2A52 AS'2 

zmoas zprcp 



^moas 'S'o(AS'^pj.j,p -I- AS'^j 



(15) 



The determination of Cin = <5'in/5'o is discussed in Sec. 
Id] Using Eq. [151 we obtain the metrological squeezing 
parameter Cm directly in terms of measured variances. 
Outside this section (|B|, we do not distinguish between 
C and Cmcas ■ 

The interesting fact that our evaluation of the metro- 
logical squeezing parameter Cm, proportional to the 
squared noise-to-signal ratio, is not affected by photon 
scattering into free space can be understood as follows: A 
longer and longer measurement of projection noise fluc- 
tuations produced by photon scattering into free space 
during the measurement, as described by Cp, results in a 
reduction of the measured spin noise. This would lead 
one to underestimate the spin noise by a factor (1 — Cp)^, 
Eq. [141 However, free-space scattering also reduces the 
observed contrast, such that the signal-to-noise ratio and 
the metrological squeezing parameter Cm do not depend 
on Cp. 

In the above derivation, we have assumed that all Ra- 
man scattering events transfer atoms between the two 
clock states. Raman scattering out of the clock states 
introduces a small correction to Eq. [T^l and a slight dis- 
crepancy between Var(Mi) and Var(M2). A full account- 
ing of the effects of different types of Raman scattering 
events changes [ASI]mi by less than 0.1 dB — i.e. less 
than the least significant figure in our results — at the 
optimum photon numbers for conditional spin noise re- 
duction and spin squeezing. 



C. Spin Noise in Rotated State 

In the inset to Fig. 1, we approximate the variance 
ASl^ of Sz in a state that has been rotated by an angle 
a about (S) between the squeezing measurement Mi and 
readout measurement M2 by 



ASi^ « Var(A/i 



Mo 



as: 



(16) 



where AS'f^^^^ = Var(Mi - M2)U=o/2. This is a good 
approximation for all a in the limit where ASz^^^ ^ 



AS'^f.gg and A^n 



> A5n 
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N 




100 120 140 160 



Duration of Microwave Pulse (^s) 

FIG. A6: Measurement of clock contrast. The open squares 
correspond to the initial contrast with only lock light, the 
other curves to probe photon numbers between p = 10^ and 
p = 9 X 10^ 



D. Contrast 

To verify and quantify spin squeezing (Fig. 3), we mea- 
sure the contrast of a Rabi oscillation after the appli- 
cation of probe light (Fig. IA6I) using an atom number 
A^Q = 4.0(1) X lO^^. The Rabi oscillation is driven by a 
microwave pulse of variable duration between the squeez- 
ing and readout measurements, during which time both 
the probe light and the resonator locking light are off. 

We observe a contrast loss that is linear in probe pho- 
ton number p, as well as a process that imparts shot- 
to-shot phase fluctuations (via imbalances in the intra- 
cavity probe power between the two spin echo pulses) 
and yields a reduction in |(S)| that is quadratic in p. 
We therefore fit to the data in Fig. 3 the expression 
C = Coexp{-ap - /3p^/2), obtaining a = 7(1) x lO"'^, 
/3 ^ 9(4) X 10~l^ and Co = 0.69(1). 

As discussed in Sec. |Bl we define Cin to be the contrast 
in the ensemble without squeezing (p—O) as observed in 
an ideal readout which flips no spins. By comparing the 
two probe pulses constituting our real readout measure- 
ment, which uses a total of 5 x lO'^ probe photons, we 
determine that the reduction in observed contrast due to 
photon scattering or imperfect microwave rotations dur- 
ing the readout is at most 4(2)%. Correcting for this 
effect, we obtain Cm = 0.71(2). 

The contrast measurement is performed at lower atom 
number than the noise measurements in Fig. 3 for two 
reasons. First, at large atom number, atom projec- 
tion noise augments the imbalances in intracavity probe 
power between the spin echo pulses and thus augments 
the resulting phase fluctuations; this effect can in prin- 
ciple be compensated using the result of the squeezing 
measurement, but we have not yet done so. Second, at 
lower atom number, the entire Rabi oscillation curve is in 
the linear regime of the Lorentzian resonator transmis- 



sion profile. We verify that the contrast is independent 
of atom number by also measuring at Aq = 3.5(3) x 10"* 
the portions of each Rabi oscillation curve that lie in 
the linear regime of the Lorentzian. A fit to these por- 
tions of the curve taken with p — (i.e., with lock light 
only) yields Cin ~ 0.63(6). While this measurement in- 
curs greater uncertainty than that at small atom number, 
it confirms that the contrast is invariant across an order 
of magnitude in atom number. 

VIII. FUNDAMENTAL LIMITS 

We outline here a derivation of the fundamental limit 
on squeezing in our system [lol . 11 1. A more complete 
treatment is given by Madsen and M0lmer [l3| . 

An ideal, Heisenberg-area-preserving measurement 
which adds no noise (i.e. a measurement where the prod- 
uct of squeezing and antisqueezing is unity) can reduce 
the variance of A^ = 2Sz by an amount inversely pro- 
portional to the broadening it imparts in phase. For our 
dispersive optical measurement, this broadening comes 
from the photon shot noise uncertainty on the AC Stark 
shift due to the probing light. For a given transmit- 
ted probe photon number p, the maximum photon-shot- 
noise-induced phase broadening is achieved by probing 
on cavity resonance and allows the normalized 5*2 noise 
to be reduced to 



Ml 



l + Nopcj^ls' 



(17) 



where (pcs — (j^orjcs/rjQ is the effective phase shift per 
transmitted photon (see Sec. |D]). However, our real mea- 
surement adds noise because photon scattering can flip 
the atomic pseudo-spins. Let Psc denote the probability 
for an atom in a superposition of states |1) and |2) to 
scatter a photon per probe photon transmitted through 
the resonator. Then in the large- detuning limit 5 3> F, 

^ Afo^eV^sc = A^O Psc = 2Afo77eff, (18) 

since (2(5/r)^Psc = ^rjcs is the single-atom resonant op- 
tical depth (see Sec. E]). Thus, Eq. [17] can be expressed 
in differential form as da^ /dp = — 2Afo?7cff-Psc(o'^)^. The 
only scattering events which contribute noise are Raman 
scattering events occurring with probability PRam per 
transmitted probe photon. We assume, for simplicity, 
the worst-case scenario that each Raman-scattered pho- 
ton transfers an atom from one clock state to the other. 
Adding the noise contribution from these spin flips and 
neglecting absorption (which for our l arg est atom num- 
ber No = 33000 is 0.6%), one obtains ^ 

da'^ 

— = -2No7]esPsc{<J^)^ + 4PRam. (19) 

dp 
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For Norjeff ^ 1, the normalized spin noise cr^ has a min- 
imum of 



-n.n = j^^. (20) 

y ^^oVcs ^sc 

In our system Pgc ~ S.OPRam- Thus, for our largest atom 
number A^o = 3.3 x 10'' and our probe polarization and 
detuning, where Norjcs = 3100, the minimum achievable 
normalized spin noise is tr^jn = — 18 dB. 

So far, we have neglected the fundamental contrast loss 
due to scattering. Including this effect fl^, one finds 

2 

that cr^jjj is reached at pPRam = ^^ln(^^) = 0.012 
for our system. The associated contrast reduction is only 

1-C ^ - v/J"R.ay,lPRay,2+-PR.am) ^ 0.012, 

where Pr^v f is the Rayleigh scattering probability in 
state \F) PJ. Thus, one also obtains a metrologi- 
cal squeezing parameter Cm,min = <7min/^^ — 
Reaching this fundamental limit would require detection 
of all the information leaving the resonator (e.g. by per- 
forming a phase measurement on resonance |13| ) with 
perfect quantum efficiency and photon-shot-noise-limited 
sensitivity. 
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